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1. Introduction 

Neutron stars can be considered as natural laboratories for studying the properties of 
matter under extreme physical conditions (see, e.g., pQ). In their cores the density p 
may exceed 10^® g cm“^ and the temperature T lies typically between 10^ and 10® K, 
but their thermal emission depends on the properties of the outer envelopes with lower p 
and T. The temperature decreases in the heat-blanketing envelopes (at p < 10^ g cm“^) 
to T ~ Teff near the radiative surface (at p ~ 10“^-10® g cm“^, depending on the stellar 
parameters), where Tgfj is the effective surface temperature related to the thermal flux 
through the Stefan’s law. In these envelopes, the magnetic held strength B may reach 
10^® G; most often B ~ 10^^ — 10^^ G. To calculate spectra of neutron-star emission 
that can be compared with the observations, one should take into account nonideality 
of the plasma in the envelopes and the effects of strong magnetic helds on the properties 
of the plasma and electromagnetic radiation. 

In the cases where ~ a few x 10^-10® K (characteristic of middle-aged neutron 

stars, i.e., of age ~ yr) and B < 10^^ G, the spectrum of thermal emission from 
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a neutron star forms in an stellar atmosphere. In a strong magnetic field, the state of 
matter can change significantly (for review, see, e.g., [2]). For instance, at Tgfj ~ 10® K, 
an atmosphere composed of hydrogen might be treated as fully ionized if the magnetic 
field were zero, but there can be a significant fraction of bound H atoms, if -B > 10^^ G. 
At lower T (characteristic of older stars) and sufficiently high B, magnetic condensation 
may occur, resulting in formation of a solid or liquid metallic surface, composed of a 
strongly coupled Coulomb plasma and covered by a thin gaseous atmosphere. 

Although the possible importance of these effects has been realized long ago, 
they were not included in neutron-star atmosphere models until recently (e.g., [21 and 
references therein). Early considerations of partial ionization in the magnetized neutron- 
star atmospheres (e.g., HEl; also reviewed briefly in [B]) were not very reliable because 
of oversimplified treatments of atomic physics and nonideal plasma effects in strong 
magnetic fields. However, in the last decade the studies of dense nonideal plasmas in the 
strong magnetic fields made a spectacular progress, that now allows one to calculate the 
spectrum and polarization of electromagnetic radiation formed in an extended partially 
ionized hydrogen atmosphere or emitted from a condensed metallic surface of a magnetic 
neutron star. In this paper, we briefly review these achievements. 

2. Equation of state and opacities of partially ionized atmospheres with 
strong magnetic fields 

Motion of electrons and ions in a magnetic held is quantized in the Landau orbitals. In 
neutron-star atmospheres, the held is usually strongly quantizing for the electrons (see, 
e.g., companion paper [7| for basic dehnitions). When B ^ 10® G, the properties of 
atoms and molecules change drastically |2]. In particular, the center-of-mass motion of 
an atom is coupled to the internal degrees of freedom. This ehect is important for both 
the equation of state (EOS) and radiative opacities of the plasma. In this section, we 
briehy review these ehects for the case of the hydrogen plasma, the only one that has 
been studied in detail up to now. This case is quite realistic, because even a small amount 
of hydrogen on the surface of the neutron star is sufficient to hll the entire atmosphere. 
For instance, for a “canonical” neutron star (with mass M = 1.4 solar masses and radius 
/? = 10 km), the total mass of the hydrogen atmosphere is ~ lO^^pTg g, where p is the 
bottom atmosphere density in g cm“^ and Tg = Teg/lO® K. This mass of hydrogen 
can be accreted from the interstellar medium or from a nebula surrounding the star, 
or produced in spallation of heavy nuclei by highly energetic plasma discharges in the 
polar-cap regions of a pulsar. 

2.1. Ionization equilibrium and equation of state 

The EOS for partially ionized hydrogen at 10^^ G< B < 10^^ G was calculated and 
discussed in [H]; recently it was extended to the case 10^^ G< B < 10^® G [2] and 
tabulated ini CD]. The treatment is based on the free energy minimization. The 
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Figure 1. Domains of partial ionization in the p — T plane for B = 10^^, 10^^ and 
10^^ G. The contours delimit the domains where the atomic fraction xn < 0.1% (to 
the right of the dotted lines), 0.1%< xh < 1% (between the dashed and dotted lines), 
1%< xh < 10% (between the solid and dashed lines) or xr > 10% (to the left of the 
solid lines). 

free energy model is in essence a generalization of the nonmagnetic model of Saumon 
and Chabrier HacH!. We consider a plasma composed of iVp protons, W electrons, 
Nn hydrogen atoms, and iVmoi molecules in a volume V, the number densities being 
Uj = Nj/V. The Helmholtz free energy is written as 

F = XJ + + Fr + Fg + f;r. (1) 

where and Fff'^ are the free energies of ideal gases of the electrons, protons, 

and neutral species, respectively, Ff^ takes into account the Coulomb plasma nonideality, 
and Fff'^ is the nonideal contribution which arises from interactions of bound species 
with one another and with the electrons and protons. Ionization equilibrium is given by 
minimization of F with respect to particle numbers under the stoichiometric constraints, 
keeping V and the total number of protons (free and bound) constant. The formulae for 
each term in (HD are given in The employed minimization technique is similar to 

the one presented in HU for the nonmagnetic case. 

Once the free energy is obtained, its derivatives over p and T and their combinations 
provide the other thermodynamic functions. These functions are affected by the strong 
magnetic held. In the domain of strongly quantizing held (p < ps, see [Zj), the electron 
degeneracy is reduced and the EOS becomes softer than it would be at i? = 0. If the 
held is weakly quantizing (p > ps), the second-order thermodynamic functions oscillate 
with increasing p around their held-free values (for details, see [SI 1101) • 

Figure n shows the domains of partial ionization in the p — T plane for diherent 
values of B. The higher B, the greater T at which the bound species are important. 
The calculated atomic fractions xh are needed to obtain the radiative opacities in the 
atmosphere (see ^2.21 below). 

Our model becomes less reliable at relatively low T where molecular chains H„ can 
be important |2]. Generally, this occurs within the p — T domain where xh ^ 0.1 (i.e.. 
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Figure 2. Radiative opacities of the plasma versus photon energy at B = 10^^ G, 
p = 1 g cm“^ and T = 10® K, for two values of the angle 9b between the propagation 
direction and the magnetic field (left: 9b = 10°; right: 9b = 60°), for the ordinary 
(chain lines) and extraordinary (solid lines) normal mode of radiation. The case of 
partially ionized hydrogen plasma (thick lines) is compared with the approximation of 
fully ionized plasma (thin lines). 


to the left of the solid lines in Fig. 

2.2. Radiative opacities, polarization and spectra 

A magnetic held affects the properties of electromagnetic radiation in plasmas (e.g., 
nil)- At photon energies hw much higher than hcppi = ~ 28.7 eV, 

where cppi is the electron plasma frequency and p is the density in g cm“^, radiation 
propagates in the form of two so-called normal modes. These modes have different 
polarization vectors Cj and different absorption and scattering coefficients, which depend 
on the angle 6b between the radiation propagation direction and the magnetic held. The 
two modes interact with each other via scattering. 

In a partially ionized atmosphere, the opacity for each normal mode is contributed 
by electrons, ions, and bound species. The scattering cross section includes contributions 
from the electrons and protons: a® = cr®’® -|- (the subscript a = 0, ±1 denotes a 
basic photon polarization - longitudinal, right or left circular - with respect to the 
magnetic held direction). The absorption cross section includes contributions from 
absorption by plasma electrons and protons (mainly by free-free transitions due to the 
electron-proton collisions, cr®), transitions between discrete states of an atom (bound- 
bound absorption, ajY) and photoionization (bound-free, crY)- Thus, for the hydrogen 
atmosphere, we can write ~ ^H(o'Y’ + crY) + (1 — a^n) cr®. For the hydrogen in a strong 
magnetic held, the cross sections were studied in US], (ty in uni, and a® in mg. 

The bound species ahect the dielectric tensor of the medium and hence the 
polarization properties of the normal modes. These ehects were studied in ma, with 
use of the Kramers-Kronig relation between the imaginary and real parts of the plasma 
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Figure 3. Spectral flux as a function of photon energy E. Left: The case of a 
partially ionized hydrogen atmosphere (solid line) at B = 10^^ G (field normal to the 
surface) and T^s = 10® K is compared with the approximation of fully ionized plasma 
(dashes) and with the blackbody spectrum (dots). Right: The cases of condensed Fe 
(B = 10^^ G, short dashes; 10^® G, dots) and H (i? = 10^"* G, long dashes) surfaces, 
at temperature T = 10® K, compared with the blackbody spectrum (solid line). 


polarizability. The full account of the coupling of the quantum mechanical structure of 
the atoms to their center-of-mass motion across the magnetic held proved to be crucial 
for the correct evaluation of the plasma polarization properties and opacities. 

An example of the opacities is shown in Fig. El In this case, the ionization degree is 
high (xh = 0.016), but the account of the bound species is still important: in the photon 
energy ranges corresponding to the bound-bound and bound-free atomic transitions, the 
model of a fully ionized plasma underestimates the opacities by a factor of ~ 10. 

These results are used in modeling the spectra of radiation from the neutron-star 
atmospheres umiHi. An example shown on the left panel of Fig. (HI clearly demonstrates 
the importance of the partial-ionization effects. 

3. Radiation from the condensed snrface in a strong magnetic field 

The notion that an isolated magnetic neutron star has a condensed surface was hrst put 
forward by Ruderman PI. who considered the Fe surface. Lai & Salpeter [20] (see also 
|2j) studied the phase diagram of strongly magnetized hydrogen and showed that, if T^s 
is below a critical value (which is a function of B), the atmosphere can undergo a phase 
transition into a condensed state. The thermal emission from the magnetized Fe surface 
was studied in jznEaES] and most thoroughly in paper j21]. For a smooth condensed 
surface, the overall emission is reduced from the blackbody by less than a factor of 2, 
whereas for a rough surface, the reduction can be negligible (see 123 ). The spectrum 
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exhibits modest deviation from blackbody across a wide energy range, and shows mild 
absorption features associated with the ion cyclotron frequency and the electron plasma 
frequency in the condensed matter. Examples of the spectrum for different models of 
the surface are shown on the right panel of Fig. El 

4. Concluding remarks 

We have briefly reviewed the main effects of strong magnetic helds on the EOS, opacities 
and properties of electromagnetic radiation in the surface layers of neutron stars. In 
order to calculate realistic spectra of thermal radiation from neutron stars, one must 
carefully take these effects into account. We expect that comparison of the calculated 
models with observations will help to improve the constraints on neutron star parameters 
(for example, their effective temperatures, see [Zj) and thus provide a powerful tool for 
testing the theories of superdense matter in neutron star interiors. 
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